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Direct  simulations  of  homogeneous  turbulence  have,  in  recent  years,  come  into  widespread 
use  for  the  evaluation  of  models  for  the  pressure-strain  correlation  of  turbulence.  While  work 
in  this  area  has  been  beneficial,  the  increasingly  common  practice  of  testing  the  slow  and 
rapid  parts  of  these  models  separately  in  uniformly  strained  turbulent  flows  is  shown  in 
this  paper  to  be  unsound.  For  such  flows,  the  decomposition  of  models  for  the  pressure- 
strain  correlation  into  slow  and  rapid  parts  is  ambiguous.  Consequently,  when  tested  in 
this  manner,  misleading  conclusions  can  be  drawn  about  the  performance  of  pressure-strain 
models.  This  point  is  amplified  by  illustrative  calculations  of  homogeneous  shear  flow  where 
other  pitfalls  in  the  evaluation  of  models  are  also  uncovered.  More  meaningful  measures  for 
testing  the  performance  of  pressure-strain  models  in  uniformly  strained  turbulent  flows  are 
proposed  and  the  implications  for  turbulence  modeling  are  discussed. 
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tions  in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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1.  INTRODUCTION 


Direct  numerical  simulations  (DNS)  of  turbulent  flows  have  provided  a  powerful  new 
tool  for  the  testing  and  screening  of  turbulence  models.  Turbulence  quantities  that  are  not 
directly  measurable  can  now  be  computed  accurately  -  for  a  variety  of  benchmark  turbulent 
flows  -  in  far  more  detail  than  that  which  can  be  extrapolated  from  physical  experiments. 
The  pressure-strain  correlation,  which  plays  a  pivotal  role  in  the  formulation  of  Reynolds 
stress  transport  models,  represents  a  prime  example  where  DNS  data  bases  have  provided  a 
wealth  of  interesting  new  information.1  During  the  past  decade  -  beginning  perhaps,  with 
the  work  of  Rogallo2  -  models  for  the  pressure-strain  correlation  have  been  tested  using 
DNS  data  bases  for  homogeneous  turbulent  flows.  In  the  most  recent  such  work,  it  has 
become  a  common  practice  to  decompose  models  for  the  pressure  strain  correlation  into 
slow  and  rapid  parts  and  to  compare  the  model  predictions  for  these  parts  separately  with 
DNS  data  bases.3-5  This  approach  seemed  attractive  since  it  would  provide  a  gauge  on  the 
performance  of  each  part  of  a  model  for  the  pressure-strain  correlation.  Prior  to  the  advent 
of  direct  simulations  of  turbulence,  such  comparisons  were  not  possible  since  the  slow  and 
rapid  parts  of  the  pressure-strain  correlation  cannot  be  measured,  even  by  indirect  means. 
However,  in  uniformly  strained  turbulent  flows,  this  procedure  can  lead  to  misleading  results. 
The  main  purpose  of  the  present  paper  is  demonstrate  this  point  along  with  other  pitfalls  in 
the  evaluation  of  models. 

It  will  be  shown  that  the  decomposition  of  models  for  the  pressure-strain  correlation 
into  slow  and  rapid  parts  is  ambiguous  for  uniformly  strained  turbulent  flows  where  the 
mean  velocity  gradients  are  constant.  Hence,  when  the  commonly  assumed  decomposition 
is  implemented  and  comparisons  are  made  with  DNS  results,  an  overly  pessimistic  and 
misleading  assessment  of  the  performance  of  a  pressure-strain  model  can  be  arrived  at.  This 
point  is  demonstrated  theoretically  and  then  is  illustrated  computationally  by  comparisons  of 
the  predictions  of  three  recent  pressure-strain  models  (the  models  of  Shill  and  Lumley,6  Fu, 
Launder  and  Tsclepidakis7  and  Speziale,  Sarkar  and  Gatski8)  with  the  direct  simulations 
of  Rogers,  Moin  and  Reynolds9  for  homogeneous  shear  flow.  It  will  thus  be  argued  that 
in  uniformly  strained  turbulent  flows,  models  for  the  pressure-strain  correlation  should  be 
tested  as  a  whole  rather  than  tested  in  their  slow  and  rapid  parts.  More  objective  alternative 
means  for  doing  so  will  be  proposed.  Of  course,  the  slow  and  rapid  parts  of  models  for  the 
pressure-strain  correlation  can  still  be  tested  separately  in  the  limit  of  relaxational  turbulent 
flows  and  in  the  rapid  distortion  limit,  respectively.  These  issues  will  be  discussed  in  detail 
in  the  sections  to  follow. 


2.  ANALYSIS  OF  THE  PRESSURE-STRAIN  CORRELATION 


We  will  base  our  analysis  on  the  incompressible  Navier-Stokes  equations 

dm  dvi  dP  2  m 

m+v^r^,+  '  () 

^  =  o  (2) 

axi 

where  ux-  is  the  velocity  vector,  P  is  the  kinematic  pressure,  and  u  is  the  kinematic  viscosity 
of  the  fluid.  The  velocity  and  pressure  are  decomposed  into  ensemble  mean  and  fluctuating 
parts,  respectively: 

Vi  =  Vi  +  rq,  P  =  P  +  p.  (3) 

An  evolution  equation  for  the  fluctuating  velocity  u,  is  obtained  by  subtracting  the  ensemble 
mean  of  (1)  from  (1)  itself  to  yield10 

dut  _  dui  du{  dv,  dp  2  ,  d  _ N  ,<x 

W  +  v’d7, =  -u^  -  u’d7,  -  fe  +  "v  +  a7,(n‘Ul)  (4) 

which  is  solved  subject  to  the  incompressibility  constraint 


We  will  restrict  our  analysis  to  homogeneous  turbulence  where  we  have 

H  -  g^> 

TquJ  =  uiuj(t). 

Hence,  for  homogeneous  turbulence,  (4)  reduces  to  the  form 

du,  dui  dui  dp  2 

dt  +  GikXkdx]  ~  UjdXj  G,]  J  dxi  +  v  *• 

By  taking  the  divergence  of  (8),  the  pressure  can  be  solved  for,  i.e., 

r^2_  du‘  du k  9ut 

y  P  -  o  r\ 

dxk  dx\  dxk 

which,  in  an  unbounded  domain,  has  the  general  solution 


"5/7/ 


X  —  X 


nr  dul  ,  dllldui\  ,3  . 

*m*i +  ojtPj '' J  ■ 
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Consequently,  (8)  can  be  written  as  an  integro-difFerential  equation  as  follows: 
dui  du i  du,  2 

r\,  Gjk^k  a  VX 

at  ox-,  oxj 


i  rj  r  x,  -  x-  f  du ;  du;  duj  \  3 

+iijjj  {2Gk,te- + d^spj d  x  ■ 


When  done  in  Fourier  space,  the  process  of  obtaining  Eq.  (11)  is  equivalent  to  projecting 
the  pressure  out  of  the  problem.11 

The  fluctuating  pressure  is  decomposed  into  slow  and  rapid  parts,  respectively,  as 
follows12,13 

p  =  p(s)  -f  (12) 

where  the  slow  pressure  p ^  and  the  rapid  pressure  p ^  are  solutions  of  the  Poisson  equations 

yys)  =  -  dxik  du‘.  mi 


dxi  dxk 


VVfi)  =  -2  Gk,p-.  (14) 

OXk 

For  homogeneous  turbulent  flows,  this  decomposition  is  unique  (for  general  inhomogeneous 
turbulent  flows  it  is  not  unique  due  to  problems  with  the  boundary  conditions).  By  making 
use  of  (12),  the  pressure-strain  correlation  can  be  decomposed  into  slow  and  rapid  parts  in 
a  straightforward  manner: 


dui  duj 
— -  +  • — - 
dxj  dxi 


=  p(s) 


du,  duj 
dx-j  dx, 


+  p(*> 


dui  diij 
dxj  dxi 


or,  equivalently, 


n,,  =  njf  +  njf* 


where  ntJ,  denotes  the  pressure-strain  correlation.  Here,  the  slow  and  rapid  parts  of  the 
pressure-strain  can  be  written  exclusively  in  terms  of  the  fluctuating  velocity  as  follows: 

—  oo 

»r  ■ 

As  a  direct  consequence  of  (11),  it  follows  that  solutions  of  the  Navier-Stokes  equations 
for  the  fluctuating  velocity  field  u,  in  homogeneous  turbulence  are  of  the  general  form  (see 
Appendix  A) 


Ul(x,t)^^[Gk,(t');  x,<],  t'  e  (0,f) 


(19) 


where  T[  •  ]  denotes  a  functional  (i.e.,  a  quantity  determined  by  the  history  of  a  function). 
From  (19),  it  is  clear  that  any  moments  constructed  from  it,  will  be  nonlinear  functionals  of 
the  history  of  Gki(t').  More  precisely,  from  (17)-(18)  and  (19)  it  follows  that  both  11-^  and 
Iljf  can  depend  nonlinearly  on  Gki  at  retarded  times  as  follows: 


n|f  (i)  =  A,,  [Gw(i');  (1,  f'e(0,i) 

n!f(t)  -  u„u [Gu(t'y,  t}Gu(t),  i'  e  (0,0 

where  the  correlation  Mijki  is  defined  by 


1  duf 
~  xl  dxl 


dui 

dxj 


+ 


d3x*. 


(20) 

(21) 

(22) 


Here,  n,^  ’  and  II b  '  are  functions  of  time  alone  since  we  are  considering  homogeneous  tur¬ 
bulence  for  which  an  ensemble  mean  is  assumed  to  be  equivalent  to  a  spatial  average. 

In  the  standard  models  for  (20)-(21)  that  have  been  used  in  the  formulation  of  second- 
order  closures,  explicit  history  effects  have  been  neglected  and  the  time  dependence  t  has 
been  parameterized  through  uiuj(t)  and  the  dissipation  rate  e(t)  -  the  only  two  turbulence 
correlations  that  are  typically  generated  in  the  solution  of  these  models.  With  these  two 
assumptions,  dimensional  considerations  and  the  fact  that  is  traceless  then  yields  models 


of  the  form12  14 


where 


n,;  =  eAij(  b)  +  KMiM  b)|| 


I<  =  -u«u, 

dui  dui 
£  V dijdxj 


bij 


U{Uj 

Ik 


(23) 

(24) 

(25) 

(26) 


are,  respectively,  the  turbulent  kinetic  energy,  the  turbulent  dissipation  rate,  and  the 
anisotropy  tensor.  Due  to  the  way  in  which  (23)  was  postulated,  it  has  become  common  to 
assume  that  models  of  this  form  can  be  decomposed  into  slow  and  rapid  parts  as  follows: 


n|f  =  £A>(b) 


(27) 


n!f  =  (28) 

Certainly,  (27)-(28)  are  consistent  with  the  limits  of  relaxational  and  rapidly  distorted  tur¬ 
bulent  flows.  More  precisely,  in  the  limit  of  a  relaxational  flow  where  for  times  t  >  0,  we  set 
0 Vk/Oxi  =  0,  it  follows  that 


n„  =  n|f  =  eAj(b) 


(29) 


1 


and  in  the  rapid  distortion  limit  where  at  time  t  =  0  we  let  GI\ o/co  — ►  oo.  it  follows  that 
for  early  elapsed  times 

rip  =  nSf>  =  ATMpW(b)^  (30) 

OX[ 

given  that  G  =  However,  for  any  uniformly  strained  turbulent  flow  where  e)Vk/e)xi 

is  constant  for  all  times,  from  (20)  and  (21)  -  which  are  a  rigorous  consequence  of  the  Navier- 
Stokes  equations  -  we  have 

since  Gki(t')  =  &vk/dx(.  Consequently,  the  decomposition  of  the  model  (2d)  into  slow  and 
rapid  parts  is  unclear  since  both  Up*  and  Ilp^  can  depend  explicitly  on  the  mean  velocity 
gradients  dvk/Ox-i.  The  only  thing  that  we  can  say  definitively  is  that  c.4q  is  the  slow 
pressure-strain  correlation  in  the  limit  of  relaxational  flows  and  that  K ,\4ijkidT>k/0xi  is  the 
rapid  pressure-strain  in  the  rapid  distortion  limit. 

This  analysis  is  particularly  relevant  for  comparison  studies  of  existing  models  for  the 
pressure-strain  correlation.  Many  of  these  models  have  been  calibrated  in  uniformly  strained 
turbulent  flows  (i.e.,  homogeneous  shear  flow).  Since  for  such  flows,  both  the  rapid  and  slow 
parts  of  the  pressure-strain  correlation  can  depend  explicitly  on  dvkjdxi ,  the  implementation 
of  the  commonly  assumed  decomposition  (27)-(28)  to  evaluate  models  can  lead  to  misleading 
results  as  we  will  show  more  clearlv  in  the  next  section. 


3.  COMPARISONS  WITH  DIRECT  SIMULATIONS 

Three  recent  models  for  the  pressure-strain  correlation  will  be  compared  with  the  direct 
simulations  of  Rogers,  Moin  and  Reynolds9  for  homogeneous  shear  flow.  These  models  are 
as  follows: 

Shih-Lumley  Model  (SL) 

Hp  =  -tfcbtj  +  -/\5q  +  12o5/\  (blkSjk  +  bjkS,k 
)  v 

—  T^bkiSfih,^  +  -(2  —  las) A  (blk\\  jk  +  bjk\\  ik) 


+  r  A '(b,ibt„,Sjm  +  blibimSim  -  '2blkSkil>i 


-3  bki$klbjj)  +  -h(b,[b  lw  II  jm  +  bjlb/m  U  ,r 


where 


(3  =  2  +  ^  exp(-7. 77 /yffUt) {72 /y/Ret  +  80.1  ln[l  +  62.4(-//  +  2.37//)]} 


F  =  1  +  9/  /  +  271 1 1 
II  —  —  bij  b{j  ^  1 1 1  =  ”  j  bj  k  bfc 

„  4  /v  2 

fle'  =  dTT 

QS  =  —  fl  +  -Z's) 

10  V  5  ) 


a5  =  —  fl  +  -F*) 

10  V  5  / 

/%  Launder  and  Tselcpidakis  Model  (FLT) 

n,j  =  120/7 \^Fe  +  1.2  (bikbkj  -  -bkibuS^ 


+  -I\Sij  +  1.2 A  (bikSjk  +  bjkSik  —  ^ bkiSki6ij 
+°^K(b,kWjk  +  bjkWtk)  +  ^K(bikbklS3l 


+bjkbkiSu  —  2  b,kSkibij  —  3bkiSkibij ) 

+  \l<  (blkbklW jt  +  bjkbklWa)  -  ~K  [siI(btkWjk 

+bjkWtkj  +  12  [blkbkiW imbmj  +  bjkbki\V  imbmj^ 


Speziale,  Sarkar  and  Gatski  Model  (SSG) 


n ij  =  -(C\e  +  C'xP)bXJ  +  C2e  (bikbkj  -  hklbkl6i ,)  +  (C3  -  c3*//| )KSi 
+CaK  (b{kSjk  +  bjkSik  —  -bkiSki6ij'j  +  C5K(btk\Vjk  +  bjkWtk) 


where  V  =  —U{U2dv{jdxj  is  the  turbulence  production  and 


Cl  =  3.4,  c;  =  1.80,  C2  =  4.2 

(41) 

c3  =  i  c;  =  1.30,  C4  =  1.25 

0 

(12) 

6'5  =  0.40,  IIb  =  btJbtJ 

(43) 

(i 


(see  Shih  and  Lumley,6  Fu,  Launder  and  Tselepidakis7  and  Speziale,  Sarkar  and  Gatski8 
for  more  details).  In  (33)  and  (39)-(40)  StJ  and  Wi}  are,  respectively,  the  symmetric  and 
antisymmetric  parts  of  the  mean  velocity  gradient  tensor  dT\/dxj  which  are  given  by 


1  /  dvi  dvj\ 

2{d7:  +  d7J 


(44) 


M(8- S) 

The  values  of  the  pressure-strain  correlation  computed  in  the  direct  simulations  of  Rogers. 
Moin  and  Reynolds9  for  homogeneous  shear  flow  will  be  compared  with  the  predictions  of 
the  SL,  FLT  and  SSG  models.  Here  the  latter  are  obtained  by  substituting  the  DNS  values 
for  bij,dvi/dxj,  I\  and  c  into  (33)-(43)  where,  in  homogeneous  shear  flow,  the  mean  velocity 
gradient  tensor  takes  the  form 

(  °  S  °  \ 

=  0  0  0  (46) 

dXj  \  o  o  o  , 


given  that  S  is  the  shear  rate.  Comparisons  will  be  made  with  the  DNS  for  t*  =  St  >  2  which 
corresponds  to  £ot/I\o  >  1.  An  eddy  turnover  time  is  allowed  for  the  artificial  early  transient 
effects  to  die  out  (the  shear  is  turned  on  at  time  t  —  0  when  the  turbulence  consists  of  a 
random  isotropic  velocity  field  with  a  square-pulse  energy  spectrum).  Since  S  is  constant, 
this  DNS  constitutes  a  uniformly  strained  turbulent  flow  as  discussed  in  the  previous  section. 

In  Figures  l(a)-(c)  the  predictions  of  the  models  for  the  slow  pressure-strain  correlat  ion  - 
based  on  the  commonly  assumed  decomposition  (27)  -  are  compared  with  the  DNS  results  of 
Rogers  et  al.9  for  homogeneous  shear  flow  (run  C128U).  From  these  results  it  would  appear 
that  the  SL  model  is  far  superior  to  the  FLT  and  SSG  models  insofar  as  the  modeling  of  slow 
pressure-strain  correlation  is  concerned.  Similarly,  in  Figures  2(a)-(c)  the  predictions  of  the 
SL,  FLT  and  SSG  models  for  the  rapid  pressure-strain  correlation  based  on  the  decomposi¬ 
tion  (28)  are  compared  with  the  same  DNS  test  case.9  Again,  one  is  tempted  to  draw  the 
same  conclusion  about  the  superiority  of  the  SL  model  in  comparison  to  the  FLT  and  SSG 
models.  However,  when  the  model  predictions  for  the  total  pressure-strain  correlation  are 
compared  with  the  same  DNS  data,  a  rather  different  picture  emerges.  In  Figures  3(a)-(c) 
the  predictions  of  the  SL,  FLT  and  SSG  models  are  compared  with  the  DNS  results  for 
the  total  pressure-strain.  Here,  the  model  predictions  are  in  much  closer  proximity  to  one 
another  and  the  relative  superiority  of  the  performance  of  the  models  is  less  clear.  The  SL 
model  is  in  the  closest  agreement  with  the  DNS  data  for  II n;  the  SSG  model  does  the  best 
for  the  ni2  component;  and  the  results  are  mixed  for  the  II22  component  although  FLT  and 
SL  are  on  balance  better  than  SSG.  From  these  results  it  is  impossible  to  judge  which  of 


t 


the  models  has  better  predictive  capabilities.  The  examination  of  other  test  cases  from  this 
DNS  data  base  is  of  little  help.  For  example,  the  same  type  of  mixed  results  are  obtained 
from  the  C12SX  run  of  Rogers  et  al.9  (see  Figures  4(a)-(c)). 

It  is  now  clear  that  for  homogeneous  shear  flow  the  direct  comparison  of  these  independent 
models  for  the  pressure-strain  correlation  with  DNS  results  for  II, j  is  not  very  illuminating. 
Furthermore,  as  shown  by  the  analysis  in  Section  2,  since  the  decomposition  of  models  for  the 
pressure-strain  correlation  into  slow  and  rapid  parts  is  ambiguous  in  homogeneous  shear  flow, 
highly  misleading  conclusions  can  be  drawn  when  comparisons  are  made  exclusively  based 
on  the  decompositions.  A  few  further  comments  are  in  order  concerning  the  implications  of 
this  point  for  turbulence  modeling.  It  could  be  argued  that  since  many  of  the  authors  who 
have  developed  models  for  the  pressure  strain  correlation  have  utilized  the  decompositions 
(27)-(28)  in  the  derivation  of  their  models,  then  why  is  it  unsound  in  these  cases  to  compare 
the  parts  of  the  model  separately?  It  must  be  remembered  that  the  rapid  parts  of  these 
models  have  usually  not  been  calibrated  separately.  Namely,  uniform  shear  flow  results 
have  typically  been  used  to  calibrate  the  “rapid”  part  of  the  model  after  the  “slow”  part 
was  calibrated  based  on  the  return  to  isotropy  problem.  Since,  in  uniform  shear  flow,  both 
the  slow  and  rapid  parts  of  the  pressure-strain  correlation  depend  explicitly  on  the  mean 
velocity  gradients  as  shown  in  (31)-(32),  it  follows  that,  when  calibrated  in  this  fashion, 
the  model  (28)  ends  up  accounting  for  both  the  rapid  pressure-strain  correlation  and  the 
strain-dependent  part  of  the  slow  pressure-strain  correlation.  This  could  explain  why  these 
models  do  not  perform  well  in  the  rapid  distortion  limit  (see  Reynolds15). 

While  the  Shih  Lumley  model  appears  to  yield  slow  and  rapid  parts  that  compare  well 
with  DNS  data  for  homogeneous  shear  flow  as  shown  in  Figures  1-2,  this  agreement  is 
fortuitous.  For  other  homogeneous  turbulent  flows,  the  agreement  of  the  individual  parts  of 
the  model  is  not  nearly  so  favorable.  In  Figures  5(a)-(c)  the  components  of  the  slow  pressure- 
strain  correlation  predicted  by  the  SL,  FLT  and  SSG  models  are  compared  with  the  DNS 
data  of  Rogallo2  for  plane  strain.  It  is  clear  from  these  results  that  none  of  the  models 
are  able  to  yield  results  in  close  range  of  this  DNS.  Similar  negative  results  were  recently 
reported  by  Shih  et  al.3  for  the  axisymmetric  expansion  as  well  as  for  plane  strain  turbulence. 
It  is  thus  clear  that  this  commonly  adopted  practice  of  testing  the  “slow”  and  “rapid” 
parts  of  the  pressure-strain  correlation  in  homogeneously  strained  turbulent  flows  should  be 
abandoned  in  favor  of  alternative  tests  based  on  the  total  pressure-strain  correlation.  This 
point  led  the  present  authors  to  model  the  total  pressure-strain  correlation,  without  using 
the  decompositions  (27)-(2S),  when  formulating  the  SSG  model.  Such  an  approach  had  been 
proposed  previously  by  Leslie.16  It  must  be  remembered  that  only  the  total  pressure-strain 
correlation  enters  into  the  Reynolds  stress  transport  equation. 
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4.  ALTERNATIVE  TESTS  OF  PRESSURE-STRAIN  MODELS  USING  DNS 


As  shown  in  the  last  section,  the  direct  comparison  of  models  for  the  pressure-strain 
correlation  with  DNS  can  be  either  misleading  or  inconclusive.  The  only  way  to  obtain  a 
more  accurate  gauge  on  the  performance  of  pressure-strain  models  in  homogeneously  strained 
turbulent  flows  is  to  compute  these  flows  with  a  full  second-order  closure  that  incorporates 
these  models  and  then  compare  with  DNS  results.  In  this  way,  the  ability  of  these  models 
to  yield  good  predictions  for  the  time  evolution  of  the  Reynolds  stresses  as  well  as  the 
equilibrium  anisotropies  can  be  assessed.  In  the  final  analysis,  the  purpose  of  these  pressure- 
strain  models  is  to  yield  a  full  Reynolds  stress  closure;  hence,  such  models  should  be  judged 
good  or  bad  based  on  their  ability  to  yield  accurate  predictions  for  the  Reynolds  stresses. 

The  Reynolds  stress  transport  equation  takes  the  form10 

•  _ dVj  _ dVi 

utUj=  -  UjUfc— ~  +  n,j  -  e,j  (47) 

for  homogeneous  turbulence.  Consequently,  in  order  to  achieve  closure,  a  model  for  the 
turbulent  dissipation  rate  tensor  etJ  is  needed  in  addition  to  a  model  for  the  pressure-strain 
correlation  11^.  This  does  cause  some  concern  since  inaccurate  model  predictions  can  arise 
from  two  sources:  the  pressure-strain  correlation  or  the  dissipation  rate  tensor.  However,  in 
virtually  all  of  the  commonly  used  second-order  closures,  the  turbulent  dissipation  rate  £  is 
modeled  as  follows: 

_  e _ dv{  e2 

e  =  ~CeXKUiUjdrj~C*~K  (48) 

in  homogeneous  turbulence  where  e  =  \sxl  and  C€\  and  Cs2  are  either  constants  or  functions 
of  II ,  III  and  Ret.  This  model  -  which  has  its  origins  in  the  work  of  Davidov1'  and  Launder 
and  co-workers18-19  -  has  been  shown  recently  by  Speziale  and  Mac  Giolla  Mhuiris20  to 
perform  well  in  homogeneous  shear  flow.  Hence,  if  versions  of  (48)  that  are  compatible  with 
the  SL,  FLT  and  SSG  models  are  chosen,  a  fair  basis  of  comparison  is  established  between 
these  pressure  strain  models.  The  form  of  the  coefficients  chosen  by  SL,  FLT  and  SSG  are 
as  follows: 

Shift  and  Lumley  Model 

Cel  =  1.2,  Ce 2  =  T  +  0.49exp(-2.83/?et~>)[l  -  0.33/»(l  -  55//)]  (49) 


Fu,  Launder  and  Tselepidakis  Model 


Cel  =  1.45,  Ce2  =  1.90 


Speziale,  Sarkar  and  Gat  ski  Model 


Cn  =  1.44,  Ce 2  =  1.83. 


(51) 


The  effect  of  variations  of  these  coefficients  will  be  discussed  briefly  later.  Both  the  SL 
and  SSG  models  are  used  with  Kolmogorov’s  assumption  of  local  isotropy 


~  3eSij- 


(52) 


In  the  FLT  model,  a  simple  algebraic  model  is  used  to  parameterize  the  anisotropy  of  dissi¬ 
pation  as  follows: 

£,J  =  ^ £\fF8ij  +  (1  -  '/F)^TTJn]  (53) 

wherein  £  is  obtained  from  (48).  This  model  is  needed  for  consistency  with  the  limit  of 
two-component  turbulence  (in  more  recent  versions  of  the  model  of  Launder  and  co-workers, 
the  anisotropic  part  of  (53)  is  combined  with  14^). 

Now,  computations  of  the  Reynolds  stresses  predicted  by  the  models  will  be  compared 
with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear.  In  Figure  6,  the  predictions 
of  the  SL,  FLT  and  SSG  models  for  the  time  evolution  of  the  turbulent  kinetic  energy  are 
compared  with  run  C128W  of  the  DNS9  (here.  /\‘  =  A'/A'o,  V  =  St).  These  results  appear 
to  indicate  that  the  SL  model  performs  the  best,  by  far.  Sliih  et  al.3  reported  very  similar 
results  for  the  SL  model  for  run  C128YV  of  the  DNS  which  they  used  to  argue  for  the 
superiority  of  this  model.  However,  it  becomes  clear  that  these  good  results  are  fortuitous 
when  other  cases  are  considered.  In  Figures  7-8,  the  model  predictions  are  compared  with 
DNS  results  for  runs  C12SU  and  C128X  of  Rogers  et  al.9  For  these  cases,  none  of  the  models 
are  able  to  predict  the  trends  of  the  DNS.  The  reason  for  this  is  that  the  initial  conditions 
of  the  DNS  are  contaminated.  More  specifically  at  time  t*  =  0  when  the  shear  is  turned 
on,  the  flow  is  seeded  with  a  random  isotropic  velocity  field  with  a  square  pulse  spectrum 
(see  Rogallo2).  Hence  it  is  unrealistic  to  expect  one-point  turbulence  models  to  predict  the 
early  time  evolution  of  this  “pseudo-turbulence”  (the  artificial  nature  of  the  early  transient  is 
best  illustrated  by  run  C128X  shown  in  Figure  8  for  which  the  turbulence  Reynolds  number 
drops  an  order  of  magnitude  during  the  first  eddy  turnover  time).  At  least  one  eddy  turnover 
should  be  allowed  to  elapse  before  making  comparisons;  this  corresponds  to 


or  St  >  2  for  the  initial  conditions  of  these  simulations.  Hence  at  St  =  2,  the  values  of 
TTTTTy  and  £  taken  from  the  DNS  which  correspond  to  a  more  physical  turbulence  with  a 
developed  energy  spectrum  are  now  used  as  initial  conditions  for  the  models.  With  these 
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more  physical  initial  conditions,  a  different  picture  emerges.  For  the  C128W  run  shown  in 
Figure  9(a)  the  SSG  model  now  yields  the  best  predictions  in  comparison  to  the  DNS  results. 
For  the  C128U  run  shown  in  Figure  9(b),  the  FLT  model  performs  the  best;  whereas  for 
the  C128X  run  shown  in  Figure  9(c)  the  FLT  and  SSG  models  perform  comparably  well. 
However,  in  all  of  these  runs,  the  SL  model  performs  poorly.  This  is  due  to  its  underprediction 
of  the  growth  rate  -  a  result  that  will  be  shown  later.  By  basing  the  calculations  on  the 
contaminated  initial  conditions  at  St  —  0,  it  would  be  erroneously  concluded  that  the  SL 
model  performs  the  best  among  these  models  in  homogeneous  shear  flow. 

Finally,  we  will  provide  a  much  more  reliable  criterion  for  judging  the  predictive  capabil- 
ities  of  pressure-strain  models  in  homogeneous  shear  flow.  All  of  these  models  predict  that 
when  an  initially  isotropic  turbulence  is  suddenly  subjected  to  a  mean  strain  then,  during 
the  early  times, 

n,-j «  g/v'N.j.  (5i) 

D 

This  result  was  first  derived  by  Crow21;  it  is  a  rigorous  constraint  that  can  be  obtained 
by  simple  symmetry  arguments.  Eq.  (54)  guarantees  that  the  models  will  perform  well  for 
early  times.  Hence,  the  prospect  that  the  models  will  perform  well  at  later  times  is  tied 
strongly  to  their  ability  to  predict  the  equilibrium  values  (btj) ^  and  ( SK/s )oc  since  these 
are  approached  within  a  few  eddy  turnover  times  and  achieve  values  that  are  independent 
of  the  initial  conditions  -  a  feature  that  is  universal  in  homogeneous  shear  flow.20  In  Table 
1,  the  equilibrium  values  predicted  by  the  SL,  FLT  and  SSG  models  for  homogeneous  shear 
flow  are  compared  with  the  DNS  results  of  Rogers  et  al.9  averaged  over  the  six  runs  discussed 
therein.  Here, 

"P\ 


(V\  o,,  x  (SK' 

(l)oo  “  "2(6li)”  (~ 
is  the  equilibrium  value  of  the  ratio  of  production  to  dissipation  and 


^oo  —  — 2(6i2)oo  —  (57^) 


(55) 


(5G) 


is  the  equilibrium  growth  rate  (for  t"  1,  /\*,e" 


~  c 


).  It  now  becomes  clear  why  the 


FLT  and  SSG  models  perform  more  favorably  compared  to  the  DNS  results  of  Rogers  et  al.9 
Both  FLT  and  SSG  yield  growth  rates  that  are  in  the  range  of  the  DNS:  similarly  the  other 
equilibrium  values  are  close  to  the  DNS  results.  However,  the  SL  model  yields  equilibrium 
values  for  btJ  and  S  I\  /  £  which  -  with  the  exception  of  (/q  i),x,  compare  poorly  with  the 
DNS  results  as  shown  in  Table  1.  In  particular,  the  growth  rate  \tyc,  is  uuderprodicted 
by  20%:  a  result  which  explains  why  the  SL  model  predictions  are  consistently  below  tin' 
DNS  data  for  the  turbulent  kinetic  energy  shown  in  Figures  9(a)-(c)  for  three  independent 
cases.  Similar  conclusions  can  be  drawn  when  comparisons  are  made  with  the  independent 
large-eddy  simulation  of  Bardina  et  al.22  shown  in  Figure  10. 


Some  comments  are  in  order  concerning  the  effect  of  the  model  constants  in  (4S)  on 
the  performance  of  the  SL  model.  For  example,  could  the  bad  predictions  of  the  SL  model 
be  due  to  problems  with  the  modeled  dissipation  rate  transport  equation?  The  answer  to 
this  question  appears  to  be  no.  Sliih  (private  communication)  has  used  a  variety  of  values 
for  C€ i  in  the  range  of  1  <  Ce\  <  1.25  for  the  calculation  of  homogeneous  turbulent  flows 
(for  homogeneous  shear  flow,  Shih  et  al.3  used  a  value  of  Ce\  «  1.2  which  prompted  us  to 
use  the  same  value  for  the  calculations  presented  herein).  It  can  be  shown  that  C£ i  —  1 
constitutes  a  bifurcation  point  of  the  t-transport  equation  where  ( SK/e )oo  -H >  oo  and  the 
exponential  growth  of  I\  and  t  is  replaced  with  an  algebraic  growth  (see  Speziale23).  Hence, 
a  significant  reduction  of  C£ i  below  1.2  will  further  degrade  the  already  bad  prediction  of 
(5A'/t)oo  made  by  the  SL  model  (see  Table  1)  and  bring  the  model  dangerously  close  to  the 
bifurcation  point  C£\  =  1.  Larger  values  of  C£ i  alleviate  the  bad  predictions  for  (SK/c) 
and  ( T,/t)<x> •  Speziale  et  al.24  computed  homogeneous  shear  flow  for  the  SL  model  using  the 
more  traditional  value  of  Ce\  =  1.44.  This  calculation  yields  ( SI\/e ),*>  ~  7  and  (V/£)lX,  ~  1.7 
which  are  much  better;  however,  it  was  at  the  expense  of  the  growth  rate  which  drops  to 
0.09  -  a  value  which  would  even  further  lower  the  predictions  of  the  SL  model  in  comparison 
to  the  DNS  results  of  Rogers  et  al.9  shown  in  Figures  9(a)-(c).  Furthermore,  for  C£\  >  1.4 
the  maximum  anisotropy  (/qi)oo  drops  to  only  approximately  one-half  of  the  value  shown  in 
Table  1  -  a  serious  underprediction.24  Hence,  it  is  clear  that  there  is  no  value  of  C£\  which 
will  render  good  predictions  for  the  SL  model  in  homogeneous  shear  flow.  Since  C€\  and  Ce 2 
only  enter  into  the  prediction  of  the  equilibrium  values  through  the  combination20 


c, 


c2 


1 


c„-rv£/„  (57) 

(which  should  be  approximately  1.8),  only  Ct  1  needs  to  be  varied.  It  therefore  follows  that 
it  is  not  possible  for  the  SL  model  to  yield  good  results  for  homogeneous  shear  flow  when 
used  in  conjunction  with  a  modeled  dissipation  rate  transport  equation  of  the  standard  form 
(48). 

Of  course,  these  results  raise  the  question  as  to  why  the  SL  model  performs  so  poorly  yet 
appears  to  be  so  favorable  in  comparison  to  the  DNS  data  of  Rogers  et  al.9  for  the  pressure 
strain  correlation  (see  Figures  3(a)-(c)).  A  possible  explanation  arises  from  Figure  3(b)  for 
the  crucial  n)2  component  which  contributes  strongly  to  the  determination  of  the  growth 
rate.20  By  the  end  of  the  simulation,  the  SL  model  predictions  are  diverging  away  from  the 
DNS  data  whereas  the  FLT  and  SSC  model  predictions  arc  converging  toward  the  DNS  data. 


5.  CONCLUDING  REMARK 


Considerable  care  must  be  exorcised  when 


DNS  data  bases  are  used  to  evaluate  the 
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performance  of  turbulence  models.  The  recent  practice  of  comparing  separately  the  model 
predictions  for  the  “slow”  and  “rapid”  parts  of  the  pressure-strain  correlation  with  DNS 
data  bases  for  homogeneously  strained  turbulent  flows  is  fundamentally  unsound.  This  is 
due  to  the  fact  that  the  decomposition  of  models  for  the  pressure-strain  correlation  into 
slow  and  rapid  parts  is  ambiguous  in  such  flows  since  both  parts  can  depend  explicitly 
on  the  mean  velocity  gradients.  As  a  result  of  the  way  in  which  models  for  the  pressure- 
strain  correlation  are  calibrated,  the  part  of  the  model  that  depends  explicitly  on  the  mean 
velocity  gradients  is  actually  accounting  for  contributions  from  both  the  slow  and  rapid 
terms  in  these  flows.  Consequently,  when  a  model  is  decomposed  into  slow  and  rapid  parts 
along  the  traditional  lines  (27)-(28),  highly  misleading  results  can  be  obtained  when  separate 
comparisons  are  made  with  DNS  results  for  the  slow  and  rapid  pressure  strain  correlations 
drawn  from  homogenous  shear  flow. 

The  individual  slow  and  rapid  parts  of  the  commonly  used  models  can  only  be  properly 
tested  in  the  limit  of  relaxational  flows  and  in  the  rapid  distortion  limit,  respectively.  None 
of  these  models  are  consistent  with  DNS  or  RDT  in  the  rapid  distortion  limit  as  recently 
demonstrated  by  Reynolds15  -  a  deficiency  which  in  all  likelihood  arises  from  the  neglect 
of  the  nonlinear  history  effects  in  (20)-(21).  Evidence  is  rapidly  accumulating  that  models 
of  the  commonly  assumed  form  (23)  cannot  accommodate  both  the  rapid  distortion  and 
moderately  strained  flow  limits.  Other  pitfalls  can  be  encountered  when  DNS  data  bases 
are  used  to  test  pressure-strain  models  in  the  limit  of  relaxational  flows  (i.e.,  the  return  to 
isotropy  problem).  The  DNS  data  bases  of  Lee  and  Reynolds25  that  are  typically  used  are  for 
Reynolds  numbers  Re\  <  10  based  on  the  Taylor  microscale.  At  such  low  Reynolds  numbers 
there  is  virtually  no  return  to  isotropy  -  an  effect  that  can  only  be  captured  by  pressure- 
strain  models  where  the  Rotta  coefficient  goes  to  2  as  Re\  becomes  of  0(1).  Consequently,  it 
is  not  surprising  that  the  SL  model  -  which  has  just  such  a  low  Reynolds  number  correction 
-  performs  better  than  the  high  Reynolds  number  SSG  and  FLT  models  in  these  relaxational 
flows  as  shown  by  Shih  and  Lumley.5  However,  when  high  Reynolds  number  experiments  are 
considered,  a  drastically  different  picture  emerges  as  shown  recently  by  Sarkar  and  Speziale21’ 
(here,  the  SSG  model  performed  better  than  the  SL  model  in  three  out  of  four  experimental 
test  cases  of  the  return  to  isotropy  problem  at  higher  Reynolds  numbers). 

In  the  final  analysis,  simplified  models  for  the  pressure-strain  correlation  such  as  (23) 
are  primarily  useful  in  so  far  as  they  are  able  to  render  reliable  predictions  for  the  Reynolds 
stresses  within  the  framework  of  an  otherwise  internally  consistent  second-order  closure.  For 
homogeneous  turbulence  this  is  best  achieved  by  satisfying  the  Crow21  constraint  and  by 
calibrating  the  model  to  yield  good  equilibrium  values  for  the  Reynolds  stress  anisotropies 
in  benchmark  cases  like  homogeneous  shear  flow  (the  former  ensures  good  predictions  for 
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early  times  whereas  the  latter  enhances  the  prospects  for  good  predictions  at  later  times). 
The  FLT  and  SSG  models  discussed  in  this  paper  perform  much  better  than  the  SL  model 
in  homogeneous  shear  flow  since  the  former  two  models  were  partially  calibrated  based  on 
experimental  data  for  this  flow.  On  the  other  hand,  the  SL  model  was  largely  formulated 
based  on  the  satisfaction  of  the  constraint  of  realizability  in  the  limit  of  two-component  tur¬ 
bulence.  The  formulation  of  models  exclusively  by  the  satisfaction  of  an  extreme  constraint 
-  exact  as  it  may  be  -  may  not  guarantee  that  a  model  will  perform  well  in  basic  benchmark 
flows.  Deficiencies  like  those  uncovered  in  this  study  can  be  overlooked  when  models  for  the 
pressure-strain  correlation  are  compared  in  isolation  with  DNS  data  bases  without  a  direct 
examination  of  their  predictive  capabilities  within  the  framework  of  a  full  Reynolds  stress 
closure. 
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APPENDIX  A 


The  equation  of  motion  (11)  for  the  fluctuating  velocity  can  be  rewritten  in  the  alternative 
form 

dui  2  du{  du{ 

dt  „V  ut  -  ujdx _  Gjtxkdx_  GijUj 

(Al) 

[]f^%  (tgM  +  d£p-\ 

|x  — x*|3\  ox'dxl) 

for  any  homogeneous  turbulent  flow.  By  making  use  of  the  Green’s  function  for  the  diffusion 
operator,  (Al)  can  be  converted  to  the  equivalent  form 

—  no  *'  J 


. ¥ 


UX^  ^  / 

oo 

+  4tt^ ///  ^(x  ~  x\Ou>(x,>0)dV 


H(x-xf,t-t')  =  — — - - ~r  exp  [  |x  —  x,|2/4i/(<  -*')]  (A3) 

2\/7n/(2  —  <')2 

is  the  Green’s  function  which  is  non-zero  only  for  t  —  t'  >  0.  For  an  ensemble  of  initial 
conditions  u;(x, 0),  one  could  envision  solving  (A2)  iteratively.  These  solutions  would  be  of 
the  general  mathematical  form 


«i(x,  0  =  /i[G«(<');  x,  t],  t'  €  (0, 0 


where  •  ]  denotes  a  functional  in  the  history  of  the  mean  velocity  gradient  tensor  Gki- 
Consequently,  the  fluctuating  velocity  -  as  well  as  any  statistics  constructed  from  it  -  will 
depend  nonlinearly  on  the  history  of  the  mean  velocity  gradients.  On  physical  grounds,  the 
duration  of  this  history  dependence  would  be  expected  to  be  of  the  order  of  an  eddy  turnover 
time. 
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Equilibrium 

Values 


SL 

Model 


FLT 

Model 


SSG 

Mode! 


DNS 

Results 


I 


(^11  )oo 

0.202 

0.208 

0.204 

0.215 

(612)00 

-0.081 

-0.146 

-0.156 

-0.158 

(622)00 

-0.195 

-0.144 

-0.148 

-0.153 

(633)00 

-0.007 

-0.064 

-0.056 

-0.062 

(SK/eU 

21.3 

6.84 

5.98 

5.70 

(V/eU 

3.43 

2.00 

1.87 

1.80 

Aqc 

0.114 

0.146 

0.145 

0.14 

Table  1.  Comparison  of  the  equilibrium  values  predicted  by  the  SL.  FLT  and  SSC  models 
with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear  flow. 
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SHEAR  FLOW  (a) 


Figure  1.  Comparison  of  the  predictions  of  the  SL,  FLT  and  SSG  models  for  the  “slow” 
pressure-strain  correlation  with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear 
flow  (run  C128U).  (a)  II^  component,  (b)  II component,  and  (c)  II 2^  component. 
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Figure  2.  Comparison  of  the  predictions  of  the  SL,  FLT  and  SSG  models  for  the  “rapid” 
pressure-strain  correlation  with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear  flow 
(run  C128U).  (a)  II  ^  component,  (b)  II ^  component,  and  (c)  11^  component. 


DNS  (C128U) 


Figure  2.  -omparison  of  the  predictions  of  the  SL,  FLT  and  SSG  models  for  the  “rapid” 
pressure-suai.i  correlation  with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear  flow 
(run  C128U).  (a)  11^  component,  (b)  II  component,  and  (c)  II  ^  component. 
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Figure  3.  Comparison  of  the  predictions  of  the  SL,  FLT  and  SSG  models  for  the  total 
pressure-strain  correlation  with  the  DNS  results  of  Rogers  et  al.9  for  homogeneous  shear 
flow  (run  C128U).  (a)  11,,  component,  (b)  n12  component,  and  (c)  fl22  component. 


PLANE  STRAIN  (a) 


<* 

Figure  5.  Time  evolution  of  the  dimensionless  “slow”  pressure-strain  correlation  in  plane 
strain:  Comparison  of  the  model  predictions  with  the  DNS  results  of  Rogallo-  (run  B2D1). 
(a)  n(ip  component,  (b)  FI^  component  and  (c)  component. 
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DNS  (B2D1) 


Figure  5.  Time  evolution  of  the  dimensionless  “slow”  pressure-strain  correlation  in  plane 
strain:  Comparison  of  the  model  predictions  with  the  DNS  results  of  Rogallo2  (run  B‘2D1). 
(a)  nS?  component,  (b)  fl^  component  and  (c)  flj^  component. 


Figure  9.  Time  evolution  of  the  turbulent  kinetic  energy  in  homogeneous  shear  flow:  Com¬ 
parison  of  the  model  calculations,  starting  from  St  =  2,  with  the  DNS  results  of  Rogers  et 
al.9  (a)  Run  C128W,  (b)  Run  C128U  and  (c)  Run  C128X. 


Figure  9.  Time  evolution  of  the  turbulent  kinetic  energy  in  homogeneous  shear  flow:  Com¬ 
parison  of  the  model  calculations,  starting  from  St  =  2,  with  the  DNS  results  of  Rogers  et 
al.9  (a)  Run  C128W,  (b)  Run  C128U  and  (c)  Run  C128X. 
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